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Thermodynamic quantities and correlation functions (CP's) of the classical antiferromagnet on the 
kagome lattice are studied for the exactly solvable infinite-component spin- vector model, D ^ oo. 
In this limit, the critical coupling of fluctuations dies out and the critical behavior simplifies, but 
the effect of would be Goldstone modes preventing ordering at any nonzero temperature is properly 
accounted for. In contrast to conventional two-dimensional magnets with continuous symmetry show- 
ing extended short-range order at distances smaller than the correlation length, r oc exp(r*/r), 
correlations in the feapome-lattice model decay already at the scale of the lattice spacing due to the 
strong degeneracy of the ground state characterized by a macroscopic number of strongly fluctuat- 
ing local degrees of freedom. At low temperatures, spin CP's decay as (SoS, ) oc 1/r^ in the range 
oo <C r <C oc T~^^^, where ao is the lattice spacing. Analytical results for the principal ther- 
modynamic quantities in our model are in fairly good quantitative agreement with the Monte Carlo 
simulations for the classical Heisenberg model, D = 3. The neutron-scattering cross section has its 
maxima beyond the first Brillouin zone; at T ^ it becomes nonanalytic but does not diverge at 
any q. 

PACS number(s): 75.50.Ee, 75.40. Cx 



I. INTRODUCTION 

Classical antifcrromagnets on kagome and pyrochlore 
lattices built of corner-sharing triangles and tetrahedra, 
respectively, are examples of frustrated systems which 
cannot ordeij-, because of the high degeneracy of their 
ground statetl and ensuing large jfluctuations. Mpnte 
Carlo (MC) simulations for kagomea and pyrochloreB lat- 
tices with the nearest-neighbor (nn) interaction J show 
a smooth temperature dependence of the heat capacity, 
C(T), in the entire temperature range. The spin corre- 
lation functions (CF's) of both models show only weak 
short-range order at T < J and decay at distances of the 
order of several lattice spacings oq. This is in a strik- 
ing contrast to the long-range order in common three- 
dimensional magnets and the extended short-range order 
[strong correlation at the distances r<^c oc exp(T*/T) 
with T* ^ J] in common two-dimensional magnets. 

Spin-wave calculations starting, from one of the or- 
dered states of the kagome latticaJ (see, also, Ref. || for 
the quantum case) yield a twofold degenerate Goldstone 
mode, as well as a zero-energy mode for all values of the 
wave vector q in the Brillouin zone, the latter reflecting 
the instability of the ground states. Mean-field approxi- 
mation (MFA) at elevated temperatures for both kagome 
and pyrochlore latticead reflects the same behavior. The 
maximal eigenvalues of the Fourier-transformed exchange 
interaction matrix are q independent for both models 
(and twofold degenerate for pyrochlore). Thus the sys- 
tem cannot choose the ordering wave vector at the mean- 



field transition temperature, , which in fact shows 

that there is no phase transition at this temperature be- 
cause of fluctuations. This results in the smooth teHU 
perature dependence of the thermodynamic qusmtitiegl3'EI 
and in the diffuse magnetic neutron scattering.El 

The degeneracy of the ground state of these models can 
be lifted by small perturbations, such as dipole-dipole 
interactions, lattice distortions, next-nearest-neighbor 
(nnn) or long-range interactions, and quantum effects. 
This may be a reason why pyrochlore antifcrromag- 
nets usually ifri^eze into a spin-glass state with lowering 
temperaturejj'u Theoretically, the most transparent way 
to lift the degeneracy is to include nnn interactions in 
the Hamiltanian.El Experiment and MC simulations on 
pyrochloresB show ordering with an unusual critical be- 
havior {(3 w 0.18) in this case. According to the spin- 
wave results of Ref. |l^ for the kagome lattice, at low 
temperatures dipole-dipole interactions favor the planar 
q = phase which is characterized by the same ordering 
pattern in each of the elementary triangles. 

A more subtle mechanism for lifting the degeneracy 
and selection of definite ordering patterns is the nonlin- 
ear interaction of spin waves for classical systems at very 
low temperatures, typically T<0.01J. For the kagome 
lattice, nonlinear effects (thermal fluctuations) favor the 
coplanar spin conflguration with the VSx short-range 
order in the case of the Heisenberg model, D =_^as. 
was suggested by the results of MC siipjulationmO'tSl 
and high-temperature series expansions.Q Extension of 
the VS X \/3 short-range order into the true long-range 
order in the limit T ^ is, however, hampered by for- 
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mation of chiral domain walls which cost no eneHgy but 
provide a gain in entropy at low concentrations. t3 The 
configuration selection at low temperatures only occurs 
if the number of spin components D is low enough. So, 
the early MC simulations of Ref. ^ for the kagome an- 
tiferromagnet showed selection of a coplanar state for 
D = 3, but no such selection for D > A. For the py- 
rochlore lattice, early simulations showed the selection of 
the collinear spiupordering for the Heisenberg model at 
low temperaturesjj although according to the recent re- 
sults of Ref. 13 this happens only for the plane rotator 



model, D = 2, and not for higher spin dimensionalities. 
The above results are in accord with general criterion for 
selection of ordered states as a function of spin and space 
dimensions for corner sharing objects, which was formu- 
lated in Ref. n3. Quantum fluctuaticms were shown to 
stabilize the v3 x ^/3 phase for S 3> l|*i|tHt-they should 
destroy ordering for low spin values S'.EHlall3'EZl One of the 
possible mechanisms for that is tunneling of the weathesy 
wane (hard hexagon) mode in the VS x \/3 structure.ll3 

It should be stressed, however, that the subtle effects 
quoted above can be easily overwhelmed by more triv- 
ial and robust ones, and they are much easier to observe 
in simulations than in experiment. The first task of the 
theory is thus to describe the principal features of classi- 
cal spin models on frustrating lattices, as, e.g., a smooth 
variation of the thermodynamic quantities in the whole 
temperature range. The simplest approach, the MFA, is 
clearly inapplicable in this case, whereas the more power- 
ful tools of the theory of critical phenomena, such as the 
renormalization group, seem to have not been yet applied 
to these lattices. 

The "next simplest" approximation for classical spin 
systems, which follows the MFA, consists in generalizing 
the Heijseal^erg Hamiltonian for the _D-component spin 
vectors! 
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|s,| = l (1.1) 



and taking the limit D oo. In this limit the problem 
becomes exactly solvable for all lattice dimensionalities, 
d, and the partition function ofJ-be. system coincidesEll 
with that of the spherical model.tJO The D ~ oo model 
possesses, however, a number of important advantages 
with respe(jt|ta,|tiie spherical one. (i) The 1/D expansion 
is possibleX- ^E^Ea including the case of low-dimensional 
systemsE3'L3 The calculations can be done conveniently 
in the framewock p©f the diagram technique for classical 
spin systemsBB'Ej'Ey (ii) Inclusion of anisotropic terms in 
Eq. (1.1) is possible, too, which allows us to clescribe or- 
dering in Icpi dimensions, including thin filmsEil and do- 
main walls (iii) In spatially inhomogeneous cases the 
D = oo model yields physically correct results, in con- 
trast to the, spherical model failing on the global spin 
constraints (iv) Below the Curie temperature Tc or in 
a magnetic field, the D = oo model describes both trans- 
verse and longitudinal CF's (Ref. that differ from 
each other, in contrast to the single CF in the spherical 
model. 



The D = oo model properly accounts for the profound 
role played, especially in low dimensions, by the Gold- 
stone or would be Goldstone modes. At the same time, 
the less significant effects of the critical fluctuation cou- 
pling leading, e.g., to the quantitatively different nonclas- 
sical critical indices, die out in the limit D oo. Thus 
this model is a relatively simple yet a powerful tool for 
classical spin systems. It should not be mixed up with tht 
TV-flavor generalization of the quantum 5* = 1/2 modeJ^ 
in the limit N oo, including its 1/iV expansionll3'L£ 
The TV-component nonlinear cr-model (see, e.g., Refs. 
as well as Ref. |3^, and ^ for the 1/N expansion) is a 
quantum extension of Eq. (^]^) in the long-wavelength 
region at low temperatures. Effective free energies for 
t he 7 1-component order parameter appear, instead of Eq. 
(1.1), in conventional theories of critical phenomena. Us- 
ing them for the 1/n expansion (see, e.g., Ref. 40) is a 
matter of taste. While yielding the same results 
critical indices as the lattice-based 1/D expansion 
it misses the absolute values of the nonuniversal quanti- 
ties. The same comment also applies to the spatially 
inhomogeneous systems in the limit D = n = oo, such as 
semi-infinite ferromagnets (cf. Refs. ^ and ^) . 

In this article the solution for the isotropic antiferro- 
magnetic infinite-component spin-vector model on the 
kagome lattice will be given. The qualitatively simi- 
lar results for the pyrochlore lattice will be presented 
in a subsequent communication. As long as the system 
studied is homogeneous, isotropic, aod-jn zero magnetic 
field, the standard spherical modeE3a can be applied, 
too. Such an approach for nonordering frustrated three- 
dimensional systems has been advocated in Ref. We 
prefer, however, to use the more general framework. 

The rest of this article is organized as follows. In Sec. 
H the structure of the kagome lattice and its collective 
spin variables are described. In Sec. Ill the formalism 



of the D — oo model is tailored for the kagome lattice. 
The diagrams of the classical spin diagram technique 
that do not disappear in the limit 13 — > oo are summed 
up. The general analytical expressions for the thermody- 
namic functions and spin CF's for all temperatures are 
obtained. In Sec. IV the thermodynamic quantities of 
the kagome antiferromagnet (AFM) are calculated and 
compared with MC simulation results in the whole tem- 
perature range. In Sec. the real space correlation 
functions are computed. In Sec. VI the neutron scat- 
tering cross section is worked out. In Sec. VII possible 
improvements of the present approach, such as the l/D 
expansion, are discussed. 



II. LATTICE STRUCTURE AND THE 
HAMILTONIAN 

The kagome lattice shown in Fig. |l] consists of corner- 
sharing triangles. Each node of the corresponding Bra- 
vais lattice (i.e., each elementary triangle in Fig. |^) is 
numbered hy i,j = 1, . . . , iV. Each site of the elemen- 
tary triangle is labeled by the index Z = 1,2,3. It is 
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FIG. 1. Structure of the kagome lattice. The elementary 
triangles are labeled by the pairs of numbers Uu , n„ according 
to Eq. (^.1), and the sites on triangles (the sublattices) are 
labeled by / = 1, 2, 3. The configuration shown corresponds to 
the coplanar \/3 x \/3 structure cha racterized by the ordering 
wave vector q^^ given by Eq. (2.3). 



convenient to use the dimensionless units in which the 
interatomic distance equals 1 and hence the lattice pe- 
riod equals 2. The triangles numbered by z, j = 1, . . . , TV 
can be obtained from each other by the translations 



(2.1) 



where r\ is the position of a site on the lattice, n„ and 
are integers, 2u and 2v are the elementary translation 
vectors (lattice periods), and 



(1,0), 



(-l/2,V3/2). 



(2.2) 



One of the most symmetric phases of the kagome AFM 
is the so-called VS x y/S phase which is shown in Fig. ||. 
This coplanar phase can be described by the complex 
"spin" Sr = Sr + — exp(«qy3'r -I- i^o), where the or- 
dering wave vector q^g can be written in three equivalent 
forms: 




(2.3) 



where w — (—1/2, — \/3/2). In this phase spins rotate by 
—240° = 120° as r changes by the lattice period 2 in each 
of the directions u, v, and w making the angle 120° with 
each other. Another realization of the -\/3 x -\/3 phase, 
in which spins rotate by —120°, is described by q^g with 
positive sign. In addition, the \/3 x \/3 phase can be 
described by appropriate combinations of different forms 
of q^ given above. 

To facilitate the diagram summation in the next sec- 
tion, it is convenient to put the Hamiltonian (1.1) into a 
diagonal form. First, one goes to the Fourier representa- 
tion according to 



q 



(2.4) 



where the wave vector q belongs to the hexagonal Bril- 
louin zone specified by the corners (±7r/3, ±7r/\/3) and 
(±27r/3,0) (see Fig. ||). The Fourier-transformed Hamil- 
tonian reads 



ll'q 

where the interaction matrix is given by 



(2.5) 



Vq = 2 J 




a = cos(u • q) 
b = cos(v • q) 
c = cos(w • q). 



(2.6) 



At the second stage, the Hamiltonian ( |2.5| ) is finally 
diagonalized to the form 



n 



nq 



(2.7) 



where V^' — 2Jj/„(q) are the eigenvalues of the matrix 
V^' taken with the negative sign, 

(2.8) 



1^1 = 1, 1^2.3 = (±Vl + 8a6c - l)/2. 
The diagonalizing transformation has the explicit form 



(2.9) 



where the summation over the repeated indices is im- 
plied and U is the real unitary matrix, = U'^ , i.e., 
U^i^ = Uin- The columns of the matrix IJ are the three 
normalized eigenvectors J7„ = (C^in, t^2n, t^3n) of the in- 
teraction matrix V: 

Un = [ac - bvn, ah - ci/„, v^^ - a^)/\/Q^, 

Qn = {vl - + [ab - cv^f + [ac - bi^„f. (2.10) 

The eigenvector Ui corresponding to the dispersionless 
eigenvalue vi — \ can be represented in the unnormalized 
form as 



Ui = [sin(w • q), sin(v • q), sin(u • q)]. 



(2.11) 



The normalized eigenvectors satisfy the requirements of 
orthogonality and completeness, respectively. 



Uln{ci)Uvn{ci)^5lv. (2.12) 



The Fourier components of the spins and the collective 
spin variables cr^ are related by 



(2.13) 



The largest dispersi onle ss eigenvalue vi of the interac- 
tion matrix [see Eq. ( p.8[ )] manifests frustration in the 
system which precludes an extended short-range order 
even in the limit T — > 0. Independence of vi of q signals 
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FIG. 2. Reduced eigenvalues of the interaction matrix, 
i^n(q) = V'q"/(2J) of Eq. plotted over the Brillouin 



that 1/3 of all spin degrees of freedom are local and can 
rotate freely. The other two eigenvalues satisfy 

z.2(q) = l-9V2, i/3(q) = -2 + -zV2 (2.14) 

at small wave vectors, = + qy 1. The eigenvalue 
i>2 which becomes degenerate with vi in the limit q ^ 
is related, as we shall see below, to the usual would be 
Goldstone mode destroying the long-range order in low- 
dimensional magnets with a continuous symmetry. The 
eigenvalue i^^ is positioned, in the long- wavelength region, 
much lower than the first two, and it is tempting to call 
it the "optical" eigenvalue. In fact, however, the eigen- 
values of the interaction matrix are not the same as the 
normal modes of the system that appear in the dynamics. 
Whereas vi gives rise to the zero-energy spin- wave branch 
corresponding to the absence of a restoring force for small 
deviations from one of degenerate ground states, i>2 and 
1/3 hybridize to the double-degenerate Goldstone mode 
with energy cx q alSiUiall wave vectors, as in conventional 
antiferromagnets.l3'u With increasing q the eigenvalue 1/2 
decreases whereas increases; at the corners of the Bril- 
louin zone they become degenerate: 1^2 = 1^3 = ^1/2. 
The q dependences of the eigenvalues i^n over the whole 
Brillouin zone is shown in Fig. 

In contrast to the smooth behavior of , the diagonal- 
izing matrix U composed of the eigenvectors Un [see Eq. 
( ^.10 )] has much more complicated structure as function 



of q. This results in the intricate behavior of the spin 
CF's and neutron-scattering cross sections at low tem- 
peratures, which will be considered in Sec. 0. Here we 
only show the nonanalytic limiting form of U at small 
wave vectors, 



-rix 



2n. 



VSrix V2 
V2, 



(2.15) 



where n = q/g. 

In the next section the equations describing spin corre- 
lation functions of the classical kagome antiferromagnet 
in the large-I? limit will be obtained with the help of 
the classical spin diagram technique. The rea ders w ho 
are not interested in details can skip to Eq. ( 3.16 ) or 
directly to Sec. IV. 



III. CLASSICAL SPIN DIAGRAM TECHNIQUE 
AND THE LARGE-D LIMIT 



The exact equations for spin correlation functions in 
the limit D ^ 00, as well as the 1/D corrections, can 
be the most conveniently obtaine d || ^v i. th the help of the 
classical spin diagram technique.EaEZle3 A perturbative 
expansion of the thermal average of any quantity A char- 
acterizing a classical spin system (e.g., A = Szi — the z 
spin component on the lattice site i) can be obtained by 
rewriting Eq. (1.1) as Ti. = Ho + Hint, where Ho is, e.g., 
the mean-field Hamiltonian, and expanding the expres- 
sion 



(A) 



1 r ^ 

- J l[dsjAeM-m, 



- 1, (3.1) 



where (3 = 1/T, in powers of Tiint- The integration in Eq. 
( |3.lD is carried out with respect to the orientations of the 
D-dimensional unit vectors Sj on each of the lattice sites. 
Averages of various spin-vector components on various 
lattice sites with the Hamiltonian Ho can be expressed 
through spin cumulants, or semi- invariants, which will be 
considered below, in the following way: 



{Sai)o — J^a, 



(3.2) 

AaApA^, 



etc., where 5ij, Sijk, etc., are the site Kronecker symbols 
equal to 1 for all site indices coinciding with each other 
and to zero in all other case s. For the one-site averages 
{i = j = k = . . .) Eq. (3^) reduces to the well-known 
representation of moments through semi-invariants, gen- 
eralized for a multicomponent case. In the graphical lan- 
guage (see Fig. |3|) the decomposition (3.2) corresponds to 
all possible groupings of small circles (spin components) 
into oval blocks (cumulant averages). The circles coming 
from Hint (the "inner" circles) are connected pairwise by 
the wavy interaction lines representing (3Jij . In diagram 
expressions, summations over site indices i and compo- 
nent indices a of inner circles are carried out. One should 
not take into account disconnected (unlinked) diagrams 
[i.e., those containing disconnected parts with no "outer" 
circles belonging to A in Eq. (3.1)], since these diagrams 
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are compensated for by the expansion of the partition 
function Z in the denominator of Eq. (3.1). Considera- 
tion of combinatorial numbers shows that each diagram 
contains the factor 1 jn^ , where is the number of sym- 
metry group elements of a diagram [see, e.g., the factor 
1/2! in Eq. ( 3.12D below]. The symmetry operations do 
not concern outer circles, which serve as a distinguish- 
able "root" to build up more complicated (e.g., renor- 
malized) diagrams. For spatially homogeneous systems, 
it is more convenient to use the Fourier representation 
and to calculate integrals over the wave vectors in the 
Brillouin zone rather than lattice sums. As due to the 



Kronecker symbols in Eq. (3.2) lattice sums are subject 
to the constraint that the coordinates of the circles be- 
longing to the same block coincide with each other, the 
sum of wave vectors coming to or going out of any block 
along interaction lines is zero. So, for our model the pair 
cumulant average of the Fourier components defined by 
Eq. (O) reads 



(3.3) 



where 8u' is the sublattice Kr onec ker symbol. The cu- 
mulant spin averages in Eq. (3.2) can be obtained by 
differentiating the generating function A(^) over appro- 
priate components of the dimcnsionless field ^ = /3HE3 



A, 



A(0=lnZo(0, 
where ^ = |||, 

Zo(0 = const X r'''/'"'^/D/2^i(0 



(3.4) 



(3.5) 



is the partition function of a Z3-component classical spin, 
and /i/(0 is the modified Bessel function. For the two 
lowest-order cumulants the differentiation in Eq. (3.4) 
leads to the following expressions: 

Aa(0-S(Oea/e, 

A.,(6 = ^(5.,-^)+S'(0^. (3.6) 
where Sap is the Kronecker symbol for spin components, 

BiO = dAiO/dC = IdMO/ lD/2-iiO (3.7) 

is the Langevin function of Z?-component classical spins, 
and B'{^) = dB/d^ (see the details in Ref. |o|). If Ho = 0, 
as is the case for our model in the absence of a magnetic 
field, the pair spin cumulant in Eq. (3.6) simplifies to the 
obvious form 



A«/3(0) = A„„(0)(5„;3, A„„(0) = 1/D. 



(3.8) 



As was shown in Ref. ^ (see also Ref. ^7|) the limit 
D oo for the spin- vector model is completely described 
by the self-consistent Gaussian approximation (SCGA), 
since all diagrams not accounted for by the SCGA van- 
ish in this limit. The SCGA consists in taking into ac- 
count pair correlations of the molecular field acting on a 
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FIG. 3. Self-consistent Gaussian approximation (SCGA) 
for classical spin systems in the nonordered state, (a) the 
Dyson equation for the spin correlation function; (b) the block 
summation for the renormalized pair cumulant spin averages. 



given spin from its neighbors, which implies a Gaussian 
statistics of molecular field fluctuations. The appropriate 
diagram sequence for the nonordered state, (sz) — 0, is 
represented in Fig. ^. Its analytical form for the square 



lattice model is given in Ref. 27, In a magnetic field 



or below in the ordering models the average spin po- 
larization (sz) 7^ appears. The additional diagrams 
and corresponding analytical expressions can be found 
in Refs. ||, |o|, and ||. The SCGA equations in the spa- 
tially inhomogeneous case and their large- 13 limit have 
been derived (and applied to domain walls) in Ref. |3|. 

In all cases above, the SCGA equations have been 
written for diagonal Hamiltonians describing the sim- 
plest one-sublattice magnets. For nondiagonal Hamilto- 
nians, such as Eq. (2.5), the matrix interaction lines, here 



— /3V^' , complicate the formalism. Simplification can be 
achieved by using the diagonalized Hamiltonian, for our 
model, Eq. (2/7). For the latter, however, the a coun- 
terparts of the one-site cumulant averages do not have a 
transparent meaning anymore since cr is a combination 
of spins on different sites and sublattices. Thus the a cu- 
mulants should be specially worked out as follows. The 
pair a cumulant (which is in our model explicitly diago- 
nal in the spin-component indices a, (3) can be rewritten 
in terms of the initial spin variables as 



/ n n \ 
\"Qa°^Qa'/0,cum 



= t//„(q)C/r„.(q')(s 



/ O.cum- 



(3.9) 



With the use of Eq. ( |3.3[ ) and the first of the relations 
(^.12D , this can be simplified to the final form 



(*^aq '-''aq' ) 0,cum. — -^aa-^^q',— q^^nn' ■ 



(3.10) 



Now, with the help of the results just obtained, the sec- 
ond diagram in the sum in Fig. ^ can be written in the 
analytical form 



A2 = C/;„(q)C//„'(q)A 



(3.11) 



where the summation over the spin-component index /3 
is implied and the quantity L/jj in the lowest (the second) 
order of the perturbation theory is given by 



I'p ni 
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(3.12) 



Here /3 — 1/T in front of V cannot be confused with 
the spin component index (3. This expression can be 
simpUfied in two ways. First, one can perform the sum 
over the index I' and use the first of Eqs. (2.12), which 
leads to 



L 



(2) 
fJl 



2! 



Tlip 



(3.13) 



Second, inverting the transformation (2.9) one can write 



(3.14) 



I'p 



Taki ng i nto account the explicit form of given by 
Eq. (2^), one can see that afte r the integration over the 
wave vector p expression (3.14) becomes independent of 
the sublattice index/. After this observation one can 
symmetrize Eq. (3.13) with respect to I. This leads to 



the vanishing of the diagonalization matrix by virtue of 
the first of Eqs. ( 2. IS ) and to the appear ance o f the factor 
1/3. Now the summation over I in Eq. (3.11) simplifies. 



and the diagonalization matrices convert, again, to S„7i'- 
The result in the second order of the perturbation theory 
has the form 



(2) 

^2 = ^aaPflL^ , 



J (2) _ ^Pf3 1 \^(ai>m^ 



3-2!7V^' 



(3.15) 



"IP 



{5nn' has been omitted) and it is independent of the 
eigenvalue index n and of the wave vector q. 

The mechanism of the simplification of diagram ex- 
pressions demonstrated above can be shown to work for 
whatever complicated diagrams. In all cases oval blocks 
represent cumulant spin averages Aqjqj . q, , as in the 
original, nondiagonalized, version of the classical spin di- 
agram technique. In all the elements connected to a given 
block summation over the eigenvalue indices n is carried 
out. The diagonalizing matrix U disappears completely 
if correlation fmictions for the a variables, 

^"(q) = §Kq<-q> = ^K'^-q): (3-16) 



N 



are considered. After the calculation of the latter the 
true spin CF's can be found from the formula 



s"'(q) = c/i4q)t^i'n(q)'^"(q) 



(3.17) 



following from Eq. ( |2.13| ). Note that (T"(q) are eigen- 
values of the correlation matrix (q) and they describe 
independent linear responses to appropriate w ave-v ector- 
dependcnt fields. As can be seen from Eq. ( 3.17 ), the 
eigenvectors describing the "normal modes" of the sus- 
ceptibility are those of the interaction matrix V^' in Eq. 



(2.5). 



The analytical expression for the a CF in the SCGA, 
which satisfies the Dyson equation shown in Fig. ||a, has 
the Ornstein-Zernike form 



1 - A„„/3V" 



(3.18) 



This expression differSj-from that obtained by Reimers 
on the mean-field basiaS by the replacement of the bare 
cumulant A^,-, = 1/D by its renormalized value A^q, de- 
termined by the diagram series Fig. The summation 
of these diagrams is documented in the most detailed way 
by Eqs. (3.16)^(3.19) of Ref. |o[ The result for A„a is 
given by the second line of Eq. ( |3.6D averaged over the 
Gaussian fluctuations of all components of the molecular 
field ^ with the dispersion defined by the quantity La ■ In 
our model, fluctuations of different components of ^ are 
independent from each other and of the same dispersion. 
La — L. Thus the quantity A^^ is diagonal and inde- 
pendent of a. In the large-D limit the multiple Gaussian 
integral determining Kaa is dominated by the stationary 
point and the result simplifles tcEII 



Aaa 



1 



D 1 + yr+si/z? 



(3.19) 



Here the dispersion L corresponding to the diagram series 
in Fig. ^ is given by the formula 



L = 



3 • 2! -^"^ 



{2t:Y 1 - kaaPV^ 



(3.20) 



generalizing Eq. (3.15). Here, the summation 
(1/A^) ... is replaced by the integration over the Bril- 
louin zone, vq is the unit-cell volume, and d is the spatial 
dimensionality. For the kagonie lattice we have vq — 2 V3 
and d = 2. The expression for L can be simplifled to 



L = 



1 

2A 



(3.21) 



where P„ is the lattice Green function associated with 
the eigenvalue n: 



Pn = Wo 



d(\ 



1 



{2TrY 1 - lacPV^ 



(3.22) 



Now one can eliminate L from Eqs. (3.19) and (3.21) 
which yields the basic equation of the large-D model. 



(3.23) 

This nonlinear equation determining Kaa as a func- 
tion pipJ-cmpprature differs from those considered 
earlieioO'EStj by a more complicated form of P reflect- 
ing the lattice structure. The form of this equation is sim- 
ilar to |tiipi| appearing in the theory of the usual spherical 
model.E3n3 The meanings of both equations are, however, 
different. Whereas in the standard spherical model a 
similar equation account for the pretty unphysical global 
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spin constraint, Eq. ( 3.23 ) here is, in fact, the normal- 
ization condition (s^) = 1 f or t he spin vectors on each of 
the lattice sites r [see Eq. (3.1)]. Indeed, calculating the 
spin autocorrelation function in the f orm sy mme trized 
over sublattices with the help of Eqs. ( 3.17 ), ( |2.12| ), and 
(3.18), one obtains 



That is, the spin- normalization condition i s au tomati- 
cally satisfied in our theory by virtue of Eq. ( 3.23 ). After 
Aaa has been found fro m th is equation, the spin CF's are 
readily given by Eqs. ( |3.1? 



and (3.17) 



To avoid possible confusion, we should menti on that in 
the paper of Reimers, Ref. ||, where Eq. ( 3.18| ) with the 
bare cumulant Aaa — has been obtained, the theo- 
retical approach has been called the "Gaussian approx- 
imation (GA)" . This term taken from the conventional 
theory of phase transitions based on the Landau free- 
energy functional implies that the Gaussian fluctuations 
of the order parameter are considered. In the microscopic 
language, this merely means calculating correlation func- 
tions of fluctuating spins after applying the MFA. Such 
an approach is known to be inconsistent, since corre- 
lations are taken into account after they had been ne- 
glected. As a result, for the kagome lattice one obtains a 
phase transition at the temperature Tc = T^^^ = 2J/D 
but immediately finds that the approach breaks down 
below Tc because of the infinitely strong fluctuations. In 
contrast to this MFA-based approach, the self-consistent 
Gaussian approximation used here allows, additionally, 
to the Gaussian fluctuations of the molecular field, which 
renormalize Aaa and lead to the absence of a phase tran- 
sition for this class of systems. The SCGA is, in a sense, a 
"double-Gaussian" approximation: The diagram series in 
Fig. jsji allows for the Gaussian fluctuations of the order 
parameter, whereas that in Fig. ||& describes Gaussian 
fluctuations of the molecular field. 

To close this section, let us work out the expressions 
for the energy and the susceptibility of the kagome anti- 
ferromagnet. For the energy of the whole system, using 
Eqs. (2.7) and ( p. 16 ), as well as the equivalence of all 
spin components, one obtains 



(H) 



N ■ 

y 



(3.25) 



To obtain the energy per spin U, one should divide this 
expression by 3A^. With the use of Eq. ( ^.18 ), the latter 
can be expre ssed through the lattice Gre en's h mction P 
of Eq. ( ^.2lD ; then with the help of Eq. ( |3.23| ) it can be 
put into the final form 



U 



2 V Aa, 



(3.26) 



The susceptibility per spin symmetrized over sublattices 
can be expressed through the spin CF's as 



1 



Xq = 



3DT 



(3.27) 



With the use of Eq. ( 3.17 ) this can be rewritten in the 
form 



1 



3DT 



:i^„2(q)a"(q), W^iq) 



I 



Ui„{q), 



(3.28) 



where the diago nalized CF's are given by Eq. ( |3.18 ). 
From Eq. ( ^.15 ) it follows that in the limit q ^ one 
has VFi = W2 = and W3 = VS. Thus the homo gcncous 
susceptibility x = Xo simplifies to 



X 



(3.29) 



As we shall see in the next section, disappearance of the 
terms with n = 1 and 2 from this formula ensures the 
nondivergence of the homogeneous susceptibility of the 
kagome antiferromagnet in the limit T — > 0. The situa- 
tion for q is much more intricate and it will be con- 
sidered below in relation to the neutron-scattering cross 
section. 



IV. THERMODYNAMICS OF THE KAGOME 
ANTIFERROMAGNET 

To put the results obtained above into the form explic- 
itly well behaved in the large- Z? limit and allowing a di- 
rect comparison with the results obtained by other meth- 
ods for systems with finite values of Z?, it is convenient to 
use the mean-field transition temperature T^^-^ — 2J/D 
as the energy scale. With this choice, one can introduce 
the reduced temperature and the so-called gap param- 
eter G according to 



T 



In these terms, Eq. ( |3.23| ) rewrites as 
eGP{G) ^ 1 



(4.1) 



(4.2) 



and det ermin es G as function of 0. Here P{G) is defined 
by Eq. (|3.2lD, where 



dq 



1 



(27r)'^l-Gi.„(q)' 



Pi = 



1 



1-G" 



(4.3) 



and the reduced eig envalu es f„(q) are given by Eq. (P^). 
The a CF's of Eq. ( 3.1^ ), which are proportional to the 
integrands of P„, can be rewritten in the form 



a"(q) = 



OG 



i-Gv^MY 



(4.4) 
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Further, it is convenient to consider the reduced energy 
per spin defined by 



U = U/\Uol Uo 



-J, 



(4.5) 



where Uq is the energy per spin at zero temperature. 
With the help of Eq. 1 0^ U can be written as 



(4.6) 



The homogeneous susceptibihtyxpf Eq. ( 3.2£ ) can be 
rewritten with the help of Eq. (2.14) in the reduced form 



G 



1 + 2G 



(4.7) 



The sense of calling G the "gap parameter" is clear 
from Eq. (4.4). If G = 1, then the gap in correla- 
tion functions closes: turns to infinity, and di- 
verges at g — !■ 0. For nonordering models, it happens 
only in the limit — > 0, however. The solution of Eq. 
(4.2) satisfies G < 1 and goes to zero at high tem- 
peratures. If 6' ^ 1, the function P is dominated by 
Pi = 1/(1 — G), whereas P3 remains of order unity and P2 
diverges only logarithmically, as in usual two-dimensional 
systems: P2 = (VS/n) ln[c/(l — G)], c ~ 1. The ensuing 
asymptotic form of the gap parameter at low tempera- 
tures reads 



G = 1- - 



V3, 3c 
— In—, 

TT 



< 1. 



(4.8) 



At high temperatures, Eq. (4^) requires small values 
of G. Here, the limiting form of P can be shown to be 
P = 1 -|- G^. The corresponding asymptote of G has the 
form 



G^i 



1 



> 1. 



(4.9) 



The numerically calculated temperature dependence of 
G is shown in Fig. IJ. Note that in the MFA one has 
G ^ 1/6 which attains the value 1 at 6* = 1. 

The temperature dependence of the reduced energy of 
Eq. (4.6) is shown in Fig. |5|. I ts asymptotic forms fol- 



lowing from Eqs. (O) and (O) are given by 



U' 



-1/9, 

-1 + {2/3)9, 



> 1 
< 1. 



(4.10) 



This implies the reduced heat capacity G = dU/d9 is 
equal to 2/3 at low temperatures, in contrast to G = 1 
for usual two-dimensional lattices in the same approxi- 
mation. The lat ter result is solely due to the term linear 
in 9 in Eq. (4^), whereas G only exponentially deviates 
from 1 at low temperatures. For the A;a(;ome lattice, there 
is a linear in 9 contribution to the gap parameter G of 
Eq. (|48|), which leads to G = 2/3. This reflects the fact 
that one of three modes in the kagome lattice [see Eq. 
(2.7)] is dispersionless, and hence 1/3 of all spin degrees 
of freedom in the system are local and free, making no 
contribution to the heat capacity. 
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e = r/rMFA 

FIG. 4. Temperature dependence of the gap parameter G 
for the kagome antiferromagnet. 
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FIG. 5. Temperature dependence of the reduced energy of 
the kagome antiferromagnet. 



The reduced variables introduced at the beginning of 
this section are very convenient for the comparison of 
the results for D = 00 with those for finite values of D, 
which are obtained by other methods. The expected dis- 
crepancies are of order 1/D which is not too much for 
D = 3. (Note that the D = ooj-aparoximation can be 
improved by the 1/D expansion.EZlc3) To compare with 
the MC simulation data of Ref. || for the heat capacity 
of the Heisenberg model we will use, instead of G, the 
true heat capacity G = dU/dT = {D/2)C [see Eqs. ( |4T| ) 
and (4^)], which in our approach tends to D/3 1 
at low temperatures. The fairlv good agreement on the 
high-temperature side of Fig. |g is not surprising, since a 
nontrivial dependence on D appears only at order 1/T^ 
for the nn correlation function and hence for the energy, 
and thus at order 1/T^ for the heat capacity [see the 
combination n + 2 = £> + 2 in Eq. (3.10a) of Ref. |. 
The reasonable agreement with the MC results at low 
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C = clUlclT 




FIG. 6. Temperature dependence of the heat capacity of 
the kagome antiferromagnet. The MC results of Ref. H for 
the Heisenberg model {D = 3) are represented by circles. 



temperatures is better than expected and can be inter- 
preted as a compensation of errors. Indeed, for finite 
values of D one should take into account the 1/D cor- 
rections to the present D = oo results. For conventional 
magnets, this leads in the first order in 1/D to the re- 
placempst-pf C ^ D/2hy C = (D - l)/2 in the limit 
T O.cJ'ES This result is exact and physically transpar- 
ent as following from the constraint |sr| = 1 in counting 
of the spin degrees of freedom; it does not change further 
at higher orders of the 1/D expansion. For the kagome 
lattice, the same counting argument suggests to replace 
D by D - 1, which would yield C ^ {D - l)/3 2/3 
for T — > 0. On the other hand, inclusion of the 1/D 
corrections reduces the degeneracy of the ground state, 
and the heat capacity should increase again. This de- 
generacy reduction manifests itself by the appearance of 
the q dependence of the correlation function a^^^ of Eq. 
(4.4). On the high-temperature side, the degeneracy of 
the largest eigenvalue of the susceptibility matrix is re- 
moved at order 1/T^ [see Eqs. (3.29) and (3.31) of Ref. 
^; the effect vanishes, however, for D ^ oo]. At low 
temperatures, the resulting heat capacity becomes 11/12 
(Ref. H), which is not far away from our result C ^ 1. 

The red uces uniform susceptibility x calculated from 
Eq. (4.7) is shown in Fig. ^ Again, our results 
are in a fairly good agreement with the MC data of 
Ref. H which are, in turn, in accord with the high- 
temperature series expansion (HTSE) results of Ref. ^ 
(not shown). Here, in contrast to the heat capacity, our 
result X — 1/3, or x = 1/(6J), at T = is exact. This 
follows from the fact that the zero-temperature suscepti- 
bilities of the classical kagome antiferromagnet have the 
same value x = 1/(6^) for all directions of the field with 
respect to the three spins on a triangle being mutually 
oriented at 120°. On the contrary, for conventional low- 
dimensional antiferromagnets, which show two-sublattice 
short-range correlations, there are D — 1 susceptibilities 
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e ^ r/r^FA 

FIG. 7. Temperature dependence of the reduced uniform 
susceptibility of the kagome antiferromagnet. The MC results 
of Ref. for the Heisenberg model {D = 3) are represented 
by circles. 



transverse to the local orientation of spins, x± — l/(2Jo)j 
and one longitudinal susceptibility X|| which vanishes in 
the zero-temperature limit. After the averaging over 
all orientations of spins one obtains the exact result 
X = (1 - l/£')/(2Jo) at T = 0, which differs significantly 
from that for the D — oo model. Taking into account 
the first-order 1 /D correction leads to a rather jaccurate 
result for xC^) the whole temperature range£3 which 
has a well-known flat maximum at T < J. Returning to 
the fcagome antiferromagnet, one can state that the 1/D 
corrections to the susceptibility are smaller than in the 
conventional case. The small maximum of x in the data 
of Ref. |l^ is probably a 1/D effect arising due to the 
increase of the longitudinal susceptibility of spins with 
temperature at low temperatures, similarly to that in 
conventional low-dimensional antiferromagnets (see Ref. 
p7| for details). 



V. REAL-SPACE CORRELATION FUNCTIONS 

The long- wavelength, low-te mpe rature behavior of the 
a corr elat io n fun ctions of Eq. (4.4) is given, according to 
Eqs. dH), (|t|), and (U), by 



3, 



3^2 



(5.1) 



where the quantity = 29/3 in cr^ defines the correla- 
tion length 



1 



261 



1/2 



(5.2) 



The appearance of this length parameter implie s tha t the 
real-space spin CF's defined, according to Eqs. (3.17) and 
(Pl),by 
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Jl' 



dq 



■'"^)C/i„(q)f/r„(q)a"(q) (5.3) 



decay exponentially at large distances at nonzero temper- 
atures. In contrast to conventional lattices, divergence of 
at — > does not lead here to an extended short-range 
order, i.e., to strong correlation at distances r < ^c- The 
zero-temperature CF's are purely geometrical quantities 
which are dominated by cr^ and have the form 



4 -^-oJj^,e^-<''~^'^'UnmMci) 

^ sin(u;-q) sin(u;--q) 

sin^(u-q) + sin^(vq) + sin^(w-q) ' 

where, according to Eq. ( p. 11 ), 

Ui = W, U2 = V, U3 = U. 



(5.4) 



(5.5) 



At large distances ftj = — i"j , the small values of q are 



important in Eq. (5.4). Thus one can expand the sines to 
the lowest order and use (u-q)^ -|- (v-q)^ -I- (w-q)^ — fq^. 
After that integration can be done analytically and yields 
the asymptotic result 



w _ 2^3 (urur)4 



2(urrij)(ui/Ti 



(5.6) 



for rij ^ 1. That is, at zero temperature spin CF's 
decrease at the scale of the lattice spacing and decay ac- 
cording t o a power law 1/r^ at large distances. The form 
of Eq. (5.6) is that of the dipole-dipole interaction in 
a two-dimensional world. Here the elementary transla- 
tion vector s u ; a ssoc iated with each of three sublattices 
[see Eqs. (5.5), (2.2), and (2.3)] play the role of dipole 
moments. 

At nonzero temperatures, an additional exponential 
decay of the correlation functions appears, wh ich is gov- 
erned by the correlation length of Eq. ( pT^). For 
<C 1, the third-eigenvalue term, n = 3, in Eq. (^.3|) can 
still be neglected, and one can use the first and second 
columns of the long- wavelength form of the diagonalizing 



matrix t//„(q), Eq. (2.15). The resulting CF s" (q) of 



Eq. (3.17), which enters Eq. (pT^), has the form 



s"'(q) = 



K2(-l + 35ii0 + 2(urq)(ur-q) 



(5.7) 



Whereas the term in the numerator yields only small 
contributions cx 6* in the real-space CF's, that in the de- 
nominator results in the additional exponentially decay- 
ing factor 



Ki{Krij) 



1, 



nn <Si 1 



TTKrij/2 e '"■'^ KM > 1, 



(5.8) 



in Eq. (5.6) 



To study real-space correlation functions at distances 
of the order of the lattice spacing and to list the partic- 
ular cases of the general formula ( |3.6[ ), it is convenient 
to en ume rate CF's by the numbers and n„ defined by 
Eq. (2J), as is shown in Fig. |l|. Thus s"^ „^ is the cor- 
relation function of the / sublattice spin of the "central" 
triangle (0,0) with the I' sublattice spin of the triangle 
translated by (n„,n„). (Note that „^ 7^ ^n^,n^.: 
general.) There is a number of several useful relations 
between correlation functions. First, the sum of the CF's 



s" n o ver ? at T = is zero by virtue of Eqs. (^) and 
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*n,0 



„22 



„33 



0, 



T = 0. 



(5.9) 



Taking into account the symmetry of the lattice, one can 
put this relation into the form of the "star rule" 



'n,0 



= 0, 



Z = 1,2,3 



(5.10) 



for the sum of the correlation functions along the direc- 
tions u, V, and w [see Eqs. (2.2) and (2.3)]. The star 
rule does not hold at nonzero temperatures, which can 
be easily seen from the HTSE for the spin CF's starting 



33 

n,0- 



from l/T^" for si]o and 5^2^ and from l/T^'^+i for s: 
More detailed analysis shows that i n the low-temperature 
region the sum in Eqs. ( |5.9| ) and ( |5.10| ) is smaller than 
^r?o by a factor of order (ku)^ ln[l/(Kn)] at the distances 
Kn ^ 1. Thus the star rule can be used with a good ac- 
curacy in the whole range 6* <C 1. An additional relation 
can be found from the condition that at zero temperature 
the sum of spins in each of the triangles is zero. Thus 
one obtains, e.g., the "triangle rule" 



J3 



0, 



T = 



(5.11) 



for all I, n„, and n„, as well as similar relations. 

The most nontrivial of the relations between spin CF's 
is the "hexagon rule" 



She 



r' ^hcx 



1)^. 



0" "rehox 



(5.12) 



for the correlators between a site r and all the sites r' 
belonging to hexagons, which are taken with alternating 
signs. If the site r itself belongs to the hexagon, the right- 
hand side of Eq. (|l2|) is nonzero and the autocorrelation 
function s^r in the sum is taken with the positive sign. 
As follows from Eq. (4.4) and the temperature depen- 
dence of the gap parameter G, the quantity Shex changes 
in this case from 1 at high temperatures to 3 at low tem- 
peratures. This very deep relation has been derived in 
Ref. 1^ from the condition that the largest eigenvalue of 
the correlation matrix ct^ is independent of q. For mod- 
els w ith fi nite D this condition and hence the hcxagoa 
rule ( |5.12D is violated only at order l/T^ of the HTSE.^ 
For our D — 00 model, given by Eq. (3.18) remains 
dispersionless at all temperatures, and the hexagon rule 
is always satisfied. 

At long distances, the zero-temperature sublattice- 
diagonal CF's in the horizontal direction, which follow 
from Eq. (p^), have the form 
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,11 



„22 



„33 



2^3 



(5.13) 



with r = 2n. One can see that relation (5.9) is satisfied. 
For the spin correlators between the first and second sub- 
lattices along the horizontal line one obtains 



„12 



2V3 



2n 
2n 



(5.14) 



Comparing Eqs. ( ^.13| ) and ( ^.14 ), one concludes that 
correlations in the D = oo model have nothing in com- 
mon with the VS X VS structure which is selected by ther- 
mal fluctuations in the Heisenberg modelJ3'Ej£3 Apart 
from the fast decay, the sign of the correlation function 
changes with each step along the line connecting the sites, 
while the coefficient alternates by the factor 2. Such a 
behavior of the sign and coefficient cannot be described 
by any ordering wave vector. Correlators involving the 
third sublattice have the form 



s^l 



„32 
*n,0 



(5.15) 



One can see tha t the above expressions satisfy the tri- 
angle rule, Eq. ( 5.11 ). In addition, the CP's along the 
three lines going through the apexes of the David stars 
in Fig. |l| read 
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„22 _ 33 



2^3 



(5.16) 
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FIG. 8. Real-space correl atio n functions s^^^ and [sl^^] at 
T = calculated from Eq. (5^). The distance unit is th e in- 
teratomic spacing. The asymptotes given by Eqs. (5TJ) and 
( [5.I4I) are shown by the dashed lines. Crosses are the results 
of Ref. |l^ for the quantum system with 5=1/2 multiplied 
by 4. 



and from other relations one obtains the CF's on the 
remote side of this hexagon: 



""1 = 3a- 1/2 = -0.038 



'2,1 



= -6a -Fl = 0.076. 



(5.20) 



To calculate the sh ort-r ange cor relation functions, one 
should use in Eqs. (5.3) or ( |3.4[ ) the full form of the 
diagonalizing matrix [/;„( q) [se e Eq. ( 2.10| )] instead of 
its long-wavelength form ( |2.15| ) and integrate over the 
whole Brillouin zone. This seems to be impossible to do 
analytically, but at T = one can express numerous CF's 
through some "fundamental" one with the help of the 
relations discussed above. So, in addition to the trivial 
results Sqq = 1, Sq\ = s^^, — ^1/2, etc., one obtains 
numerically 



'"0 



„22 
1,0 



0.1540. 



(5.17) 



After that using the star and triangle rules leads to the 
results 



„33 
^1,0 

ol3 



= -2a = -0.308 



''1,1 

''1,0 
„11 



,13 



'1,-1 
'1,-1 



1/2 = 0.346 
3a - 1/2 = -0.038 
-4a + 1/2 = 0.116 
-6a + 1 = 0.076 
3a- 1/2 = -0.038. 



(5.18) 



Now from the hexagon rule for the hexagon marked by 
the star in Fig. 



She 



'1,0 



,13 
'1,1 



'1,1 



'2,1 



'2,1 



^1,0 ~ 



(5.19) 



After that the star and triangle rules yield 

= 10a- 3/2 ^ 0.040 
sl^g ^ -36a + 11/2 -0.044 
= 26a- 4 = 0.004 
= -29a + 9/2 = 0.034. 



'2,0 



,11 



(5.21) 



This routine cannot be continued without numerically 
calculating the next fundamental CF s^^q — 0.0164. This 
would make little sense, however, because at such dis- 
tances correlation functions are already well described 
by their asymptotic forms given above (see Fig. ^). 

The results obtained above for the real-space CF's can 
be compared with those for the quantum Heisenberg an- 
tiferromagnet with S = l/2.y The latter for the CF's of 
the same type, multiplied by a factor of four to normal- 
ize the autocorrelation function by one, are also shown 
in Fig. ^. In contrast to the classical Heisenberg antifer- 
romagnet on the kagome lattice selecting the \/3 x \/3 
structure, ihe. ground state of the quantum model is 
disordered,Enlj which brings it, in a sense, closer to our 
D = 00 model. In the quantum model CF's decay faster: 
Even for the nearest neighbors one has sl% = -0.2922 
instead of —0.5 because of the zero-point motion. For the 
establishing of the large-distance behavior of the CF's in 
the quantum antiferromagnet, the numerical diagonaliza- 
tion of clusters of much larger sizes (36 sites in Ref. |l^ 
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or 27 sites in Ref. |T^) is needed, which is a tremendous 
computational problem. 

The main implication of this section is that the 
spin correlation functions for the large- model on the 
kagome lattice decay at the distances of the order of the 
lattice spacing even at T — > 0, in spite of the divergence 
of the correlation length S,c- Thus T = is not a criti- 
cal point of this model, in contrast to the conventional 
low-dimensional ferro- and antiferromagnets. Correla- 
tions developing in the low-temperature range, which, 
however, become only strong between several neighbor- 
ing spins, characterize the state of our model as a spin 
liquid. 



VI. NEUTRON-SCATTERING CROSS SECTION 

The static magnetic neutron-scattering cross section is 
proportional to the static Fourier-transformed spin CF: 



da 

In 




oc 



(6.1) 



FIG. 9. Neutron-scattering cross section from the large-D 
kagome antiferromagnet at T = (cf. Fig. ^). 



where is the component of the spin perpendicular 
to the scattering wave vector q. In our model all spin 
components are equivalent and r = r' is defined in Sec. 
||. Since the overall coefficient in Eq. (^^) contains a 
magnetic form factor and is poorly known, one can use 
the most convenient form of this coefficient and define 



da 

dn 



(6.2) 



This expression differs from the wave-vector-dependent 
susceptibility of Eqs. ( |3.27 ) and (3.28) only by the ab- 
sence of DT in the denominator, while the normal-mode 
CF's cr" are given by Eq. (4. 
vector q is not confined to the Brillouin zone 



The scattering wave 
(BZ), in 



contrast to q appearing in the calculation of the ther- 
modynamic quantities and real-space CF's. For usual 
bipartite lattices, the scattering cross section with q out- 
side the BZ is the same as with q' = q — inside the 
BZ, where k is an appropriate reciprocal-lattice vector. 
That is, in this case da/dil is repeating over the set of 
extended Brillouin zones. For periodic lattices with more 
complicated structures, the neutron cross section is still a 
periodic function of q, but the period can be larger than 
one BZ. 

For the kagome lattice in the limit T — > 0, Eq. ( |6.2| ) 
is dominated by the term with a^ = 9G/{1 — G). Using 
G = 1 — /i one obtains da/dil = Wi{c{) that is tem- 
perature independent. As follows from the contour plot 
in Fig. ^, the "unit cell" for the neutron cross section 
contains four Brillouin zones: one with a very low scat- 
tering intensity, such as the first BZ in the middle, and 
three BZ's with a highly inhomogeneous scattering pat- 
tern oriented at three different angles. The neutron cross 
section is symmetric with respect to rotations by 60° de- 
grees. It reaches its maximal value for q ~ (±87r/3,0). 



etc., and vanishes along the directions qx = 0, ±a/3(7j,, in- 
cluding q = 0. The scattering pattern in Fig. || strongly 
resembles that in the appropriate plane for the classi- 
cal Heisenberg model on the pyrochlore lattice, which 
was obtained by MC simulations in Ref. In contrast, 
the perturbative calculation for the quantum AFM model 
with S" = 1/2 on the pyrochlore latticec3 shows much less 
revealed triangular shape of the maxima of the neutron 
cross section. 

Let us now analyze the n eutro n cros s sec tion in more 
detail. At T = using Eqs. ( |2.1lD and ( ^.28D one obtains 



da _ 4sin^((7a;/2)[cos((7a;/2) 
dVl ~ 



3(V3g,/2)]2 



sin^(gx) + E± sin2(g^/2 ± Viqy/2) 



(6.3) 



where ^_|_ sums terms with both signs. Particular cases 
of this formula are the following. For = one has 



da 4sin''(ga,/4) 
'dn ^ l-F2cos2(g^/2)' 



(6.4) 



This is equal to 1/3 at the corner of the first BZ, q^ = 
27r/3, to 2 at g^; = TT (the highest slope), to 3 at g^; = 47r/3 
(the maximum], and to 8/3 at qx = 27r. Near the line 
qx = 0, Eq. ( |6.3| ) simplifies to 



^-^tan^ 
dn 2 



2tt 

71 



> qx 



For small wave vectors one obtains 
1 qlill^Hf 



da 

'dn 



96 



ql 



qx,qy < 1. 



(6.5) 



(6.6) 



The most interesting form of the neutron cross section 
is realized in the vicinity of the centers of the Brillouin 
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FIG. 10. Powder average of the neutron-scattering cross 
section from the kagome antiferromagnet at T = 0. Inset: 
MC simulations of Ref. 12 for the Heisenberg model. 



zones surrounding the first BZ in Fig. In particular, 
near q = (0, 2n/\/3) Eq. (O) takes the form 



da 



Qx^Sqy < 1, 



(6.7) 



where 6qy = Qy — 27r/-\/3. This function is nonanalytic at 
qx,Sqy = 0, and its limiting value in this point depends 
on the way of approaching to. (Such a function is difficult 
to plot: What looks like narrow paths in Fig. |9| are in 
fact infinitely thin walls.) 

At T ^ one should take into account the terms with 
n = 2, 3 in Eq. (3.2). This is especially important in the 
region where the zero-temperature neutron cross section 
turns to zero or is singular. In particular, near q = the 
quantity Wi is small due to cancellation of the leading 
t erm s and is given by the right-hand side (rhs) of Eq. 
(S.6). Similar cancellation takes place in W2, and the 
result is given by the same formula with interchanged 
q^ and qy. The leading term thus becomes the non- 
cancelled one associated with the "optical" eigenvalue, 
= 3. Thi s lead s to the uniform susceptibility given 
by Eqs. ( 3.29 ) and (4.7), that tends to a constant in the 
low-temperature limit. On the contrary, Xq with q 7^ 
behaves as 1/T at low temperatures. 

In the vicinity of the singularity point q = (0, 27r/V3), 
i.e., near the center of the Brillouin zone just above 
the first (central) BZ in Fig. ||, cancellation of the 
leading terms does not take place. Here the matrix 
(jq differs from that of Eq. ( ^.15[ ) by the redefinition 
Hy = 5qy/[q1 + {SqyYY^"^ and by the change of sign of 
the third row. Here 14^2^ is given by the rhs of Eq. (6.7) 
with interchanged qx and 5qy. Now with the use of Eq. 
(^.iD one obtains the final result 



da 

dh 



3k2 



n,qx,5qy<.l. (6.8) 



It can be seen that this expression is nonanalytic only in 
the limit T — > 0, where the correlation length defined 



cross section 
(3.2) over the direc- 



by Eq. (5.2) becomes infinite 

The powder average of the neutron 
(da/dil), i.e., the average of Eq, 

tions of q, is shown in Fig. [l^. Positions of its sin- 
gularities at T = can be found from the scattering 
pattern in Fig. ^. The first of them are located at 
q = 27r/V3 ~ 3.63 (sharp maximum), 27r « 6.28 (sharp 

^■ny/TJ?, w 9.6 (smah 

10.88 (sharp maximum), 

13.08 (sharp mini- 



mimmum 



1), ^(s^o^Twvi)^ = 

cuspy shoulder), 3 X 27r/V3 w 10 
^^(471)2-^(2^^3)2 = 2^71373 



mum ), v/(57r)2 ^ (7r/V3)2 = ^-K^fWJl w 15.81 (sharp 

maximum), 5 x 27r/-\/3 ~ 18.14 (sharp maximum), 67r k, 
18.85 (cuspy shoulder), etc. With increasing of q the be- 
havior of {da / dVL) becomes more and more irregular, and 
it very slowly approaches the value 1. The latter can be 
understood since for large values of q the average over 
the directions of q should be e qual to that over q itself. 
The latter is according to Eq. ( |6.1[ ) nothing else but the 
autocorrelation function, and the result is unity for the 
normalization of da/dfl adopted in Eq. (6.2). 

At nonzero temperatures, the sharp features of 
{da/dVL) smoothen. Their lo w-t emperature forms can be 
found with the help of Eq. (|6.8[) and are given by 



^) ^2.712 + 0. 



4n/3 {5qf 



near the first maximum, 5q = q~ 27r/V3, and 
da 



(6.9) 



— ) "Si 0.0183 - 0.2385(? + -^k^ + {6qf (6.10) 

d\l / TT 



— q — 2tt. In the high- 
is dominated by the au- 



near the first minimum, 
temperature limit, Eq. (6.1) 
tocorrelation function, and the neutron cross section de- 
fined by Eq. ( |6.2[ ) is equal to 1 for all q (an absolutely 
diffuse scattering). 

The MC simulation data of Ref. |l2| for the Heisen- 
berg model at T = 0.002 J are shown in the inset to Fig. 
p^ . At such low temperatures the system shows a ten- 
dency towards selection of the VS x ^/Z phase, and the 
corresponding Bragg-reflection peaks grow. The latter, 
according Eq. (2.3), are situated at the corners of all 
Brillouin zones in the extended BZ scheme shown in Fig. 

and their positions are shown by additional tics in Fig. 
lOj . These peaks that are superimposed on the underlying 
D = 00 structure can be traced out in the inset. Note 
that there are Bragg-reflection peaks in the vicinity of 
the D — 00 peaks, and they seem to be mixed together 
in the simulations of Ref. ^ In contrast, the first two 
minima in Fig. |l^ can be found in the inset at nearly the 
same positions, although in a strongly rounded form. 
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VII. DISCUSSION 

In the main part of this article, we have presented 
in detail the exact solution for the D = oo component 
classical antiferromagnet on the kagome lattice. The so- 
lution does not show ordering at any temperature due 
to the strong degeneracy of the ground state, and the 
thermodynamic functions behave smoothly. In contrast 
to conventional two-dimensional magnets, there is no ex- 
tended short-range order at low temperatures, and T — 
is not a critical point of the system. Although the corre- 
lation length diverges as c>c T"^/^, the power-law decay 
(sQSr) cx 1/r^ of the spin correlation functions leads to 
the loss of correlations at the scale of the interatomic 
distance. The magnetic neutron-scattering cross section 
becomes nonanalytic at T ^ but does not diverge at 
any "ordering" wave vector. 

Although the model with an infinite number of spin 
components may appear very unphysical at the first 
glance, it is in fact the second that should be applied, 
after the mean-field approximation, to any classical spin 
system. It properly describes the effect of would be 
Goldstonc modes and thus it has important advantages 
against the MFA. Properly scaled physical quantities 
show a smooth dependence on D, and in typical cases 
the large-D model proves to be a reasonable approxima- 
tion to the realistic one with D = 3. So, the results 
for the heat capacity and the uniform susceptibility ob- 
tained above are in a fairly good agreement with the 
MC simulation results for the Heisenberg model in the 
whole temperature range. This implies that the 1/D cor- 
rections to the thermodynamic functions of the kagome 
AFM, which-|(f©uld be studied within the same theoretical 
framework,cJ'c3 are suppressed by some mechanism. 

The D = oo model and the 1/D expansion seem to be 
inefficient in the cases when, due to topological effects, 
the behavior of the system abruptly changes at small 
values of D. A well-known example is the Berezinskii- 
Kosterlitz-Thouless transition which takes place in two 
dimensions for D = 2. For the kagome lattice, thermal 
fluctuations favor the \/3 x \/3 phase at low tempera- 
tures for the Hedsenjaerg model, but there is no such an 
effect for D > 3.EJ't3 On the other hand, the tendency to 
the selection of the a/S x V3 phase at low temperatures 
is already seen in the high-temperature series expansion 
of Ref. ^ for any finite value of D. Exactly how this 
mechanism becomes inefficient at low temperatures for 
D > 3, could be studied with the help of the 1/D ex- 
pansion. The latter describes lifting of the degeneracy 
of the largest eigenvalue of the correlation matrix in the 
first order in 1 /D and is applicable in the whole range of 
temperatures. 

As follows from the consideration above, the Heisen- 
berg antiferromagnet on the kagome lattice is still not the 
best model to substitute it with the D — oo model. For 
the Heisenberg AFM on the pyrochlore lattice, the large- 
D approximation can be expected to work even better 
since this model is, in a sense, more disordered, and topo- 
logical effects leading here to the selection of the coplanar 



phase arise only for D = 2.113 The formalism for the py- 
rochlore lattice in zero field is essentially the same as for 
the kagome lattice, and the corresponding results will be 
presented elsewhere. 
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